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Abstract 

We study effects of electron correlation on the transport through a small in- 
teracting system connected to reservoirs using an effective Hamiltonian which 
describes the free quasi-particles of a Fermi liquid. The effective Hamiltonian 
is defined microscopically with the value of the self-energy at u = 0. Specifi- 
cally, we apply the method to a Hubbard chain of finite size N (= 1, 2, 3, . . .), 
and calculate the self-energy within the second order in U in the electron-hole 
symmetric case. When the couplings between the chain and the reservoirs 
on the left and right are small, the conductance for even iV decreases with 
increasing N showing a tendency toward a Mott-Hubbard insulator. This 
is caused by the off-diagonal element of the self-energy, and this behavior is 
qualitatively different from that in the special case examined in the previous 
work. We also study the effects of the asymmetry in the two couplings. While 
the perfect transmission due to the Kondo resonance occurs for any odd N in 
the symmetric coupling, the conductance for odd N decreases with increasing 
N in the case of the asymmetric coupling. 
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I. INTRODUCTION 



Effects of the inter-electron interaction on the transport through small interacting sys- 
tems have been a subject of current interest. For instance, the Kondo effect expected to 
be observed in a quantum has been studied intensively from both theoretical^'! and 

experimental sides.iHl Also, the Tomonaga-Luttinger behavior expected in one-dimensional 
systemdlS has been investigated in a quantum wire.0 

For studying the quantum transport of small systems theoretically, a formulation which 
is able to treat both the interaction and interference effects in a unified way is necessary. 
Meir and WingreenEl have presented one such framework using the nonequilibrium Keldysh 
formalism. Especially, Eq. (6) of Ref. [12] is a general formula for the total current, and 
is available for various systems. While the quantum transport theory for nonequilibrium 
states is important to understand the physics in small systems, there are some ambiguities 
for determining the stationary state. So, it seems to be still meaningful to describe a 
formulation based on the linear response theory. In the Kubo formalism, the conductance 
for non-interacting electrons at T = can be expressed in terms of the Green's function at 
Fermi energyJlUll This is also true for interacting electrons if the ground state is a Fermi 
liquid-SIll The purpose of this paper is to apply a quasi-particle description of a Fermi liquid 
to the conductance of a small interacting systems. At T = 0, an effective Hamiltonian for 
free quasi-particles can be defined microscopically based on the perturbation theory, which 
reproduces the value of the Green's function at u = 0. This kind of an effective Hamiltonian 
has been introduced, for instance, for the single Anderson impurity and systems with a 
translational invarianceS One of the characteristics of the small interacting system is that 
the self-energy has off-diagonal elements. This is because generally the system consists of a 
number of resonant states and has no translational invariance. 

Specifically, we apply the method to a Hubbard chain of finite size N (= 1, 2, 3, . . .) 
connected to non-interacting leads. This system can be regarded as a model for a Mott- 
Hubbard insulator of nanometer scale. Also, for small N, the system can regarded as a model 
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for a series of quantum dots, which has been examined by advanced numerical methods 
such as the numerical renormahzation group0 and quantum Monte Carlo methodsEHl 
It may also be considered as a model for a quantum wire, which has been studied with 
the bosonization approach take into account the Umklapp scattering.&^il We calculate all 
the elements of a N x N matrix self-energy within the second order in U in the electron- 
hole symmetric case. The unperturbed Hamiltonian is taken to be connected including the 
coupling between the Hubbard chain and leads in it as that in the case of the Anderson 
impurity.i§ In the previous paper ,0 we have reported the results obtained in a special case 
v L — v R — t, where v L (v R ) is the tunneling matrix element between the chain and the 
left (right) lead, and t is the hopping matrix element of the chain. In this case, one can 
obtain the self-energy analytically owing to the property that the non-interacting system 
has the translational invariance. However, this particular feature is realized only in this 
case, and some behaviors seeing in the results seem to depend on it: the reduction of the 
conductance for even N is proportional to U A for small U, and in the limit of large even N 
the conductance tends to be finite. 

In the present work, we examine the dependence of the conductance on the coupling v L 
and v R calculating the self-energy numerically. In contrast to the special case mentioned 
above, the results obtained in weak coupling cases v L , v R < t show quite different behaviors: 
the reduction of the conductance for even iV is proportional to U 2 for small U, and the 
conductance tends to zero for large even N. Physically, the difference between the weak 
coupling cases v L , v R < t and the special case v L = v R = t is caused by whether N revels 
in the chain form distinguishable resonant states or not. The effect of the interaction is 
enhanced for the weak coupling cases. It is also shown generally beyond the second-order 
perturbation that the perfect transmission occurs for odd N when the system has both the 
inversion symmetry v L = v R (= v ) and the electron-hole symmetry. This is due to Kondo 
resonance and independent of the values of v and U. We also examine the effects of the 
asymmetry in the coupling v L 7^ v R , which disturbs the perfect transmission. We note that 
preliminary results were reported in a proceedings.^ 
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In Sec. ||, we introduce the effective Hamiltonian, and give the interpretation of the dc 
conductance and total charge displacement in terms of the free quasi-particles. In Sec. [TT1]. 
we show some properties of the quasi-particles in the electron-hole symmetric case without 
specifying details of the model. In Sec. |Vj, we apply the method to the Hubbard chain 
of finite size, and present the results obtained with the second-order perturbation theory. 
Summary is given in Sec. [V|. 

II. EFFECTIVE HAMILTONIAN AND GROUND-STATE PROPERTIES 

In this section, we introduce an effective Hamiltonian for the free quasi-particles based on 
the perturbation theory in the inter-electron interaction. In this context, the dc conductance 
and total charge displacement at T = are described by the scattering coefficients of the free 
quasi-particles. In what follows, we will discuss the formulation assuming the single-mode 
leads for simplicity. The formulation can be generalized to the multi-mode leads. 

We start with a system which consists of three regions; a finite central region (C) and 
the two non-interacting leads on the left (L) and the right (R). The central region consists 
of N resonant levels, and the interaction is switched on only for the electrons in this region. 
We assume that the two leads consist of infinite degrees of freedom, respectively, and are 
connected to the central region by the mixing matrix elements v L and v R as illustrated in 
Fig. [I]. The complete Hamiltonian is given by 

H = Hl + 7~Cr + Ti-C + 1 + 7~Lrnix , (1) 

Kl = J2 £ ( "4 - ^ 5 » ) °l c i° > ( 2 ) 
^=EE(^-' 1 ^) i c i-> ( 3 ) 

ij(zR <J 

^ = ££HS-M,)4^> ( 4 ) 

ij€C a 

{i}ec **> 
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7~Cmix — V L y C L C Ocr + C Ocr C \a J ^ ( C WlCT C iVer + C !v<7 C jV+lcr j ■ (6) 

a a 

Here cj CT (c- CT ) creates (destroys) an electron with spin a at site j, and /i is the chemical 
potential, t^, tfp and tfj are the intra- region hopping matrix elements in each regions L, R, 
and C, respectively. The coupling between the central region and two leads are described 
by the mixing matrix element 7~C m i x . We assign the labels 1, 2, . . ., N to the sites in the 
central region. Specifically, we assign the labels 1 and N to the sites at the interface on 
the left and right, respectively. Correspondingly, the label and N + 1 are assigned to the 
sites at the lead-side of the interface on the left and right, respectively [see Fig. |IJ. The 
inter-electron interaction Uj 4 j 3 -j 2 j 1 is switched on in the central region, and it has the time 
reversal symmetry: U^pi is real and = ^34 ; i2 = ^i2 ; 34 = ^42 ; 3i = f^i3 ; 24- We take 

all the hopping matrix elements to be real, and will be using units K — 1 unless otherwise 
noted. 

In the limit N = 1, the model reduces to a single Anderson impurity in which the per- 
turbation expansion is valid for all values of £7.0 Our basic idea is to apply the perturbation 
theory in TCq* following Yamada and Yosida@ regarding the central region as one big impu- 
rity. Therefore, we take the unperturbed part of the Hamiltonian to be connected by 
including the mixing term H m i x m it; 

-^(o) _ <ft L _|_ -^o_ _|_ yi mix _ (7) 

Our working hypothesis is that the ground state is changed continuously against the adia- 
batic switching-on of the interaction H.™*. The perturbation expansion can be done using 
the single-particle Green's function 

G jf (i £l ) = - J^dr ( T T c Ja (r) 4,(0) ) e<« T , (8) 

where (3 = 1/T, si = (21 + l)n/f3, c Jcr (r) = e TH Cj a e~ TH , and (• • •) denotes the thermal 
average Tr [ e _/3W ■ • • ] /Tr e~^ H . The spin index has been omitted from the left-hand side of 
Eq. (|^) assuming the expectation value to be independent of whether spin is up or down. 
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Since the interaction is switched on only for the electrons in the central region, the Dyson 
equation is written as 

G tJ (z) = Gf{z) + Gfiz) E H ,(z) G Vj {z) . (9) 
WeC 

Here G^ \z) is the unperturbed Green's function corresponding to . The summations 
with respect to I and V run over the sites in the central region, and £^(2) is the self-energy 
correction due to Ti.™ 1 . Note that G{j(z) = Gji(z) and ^(z) = £^(2) because of the time 
reversal symmetry of TC. In what follows, we will treat z as a complex variable, i.e., Gij(z) is 
the analytic continuation of Gij{iej). When the perturbation expansion is valid, the single- 
particle excitation at the Fermi energy z = i0 + does not decay at T = owing to a property 
of the Fermi liquid^ 

Im£+(0) = 0. (10) 

Here the superscript (+) is the label for a retarded function: we will use a notation Ti^{u) = 
£ij(a; ± i0 + ) H Due to this property, the Dyson equation Eq. @ for u = can be mapped 
on to a scattering problem of the free quasi-particles described by the effective Hamiltonian 

7~L qp = + Hr + He + 7~Cmix , (11) 
-He = EE( 4 S-^) C U>> (12) 

ij&C a 

-t% = -t% + ReS+(0). (13) 

The value of the Green's function corresponding to 7i qp and the one corresponding 7i are 
the same at T = 0, uj = 0. 

We now consider the conductance with the Kubo formalism. If the ground state is a 
Fermi liquid, the contributions of the vertex corrections for the dc conductance vanish at 
T = O.Ellii Therefore, the dc conductance at T = is written in terms of the Green's 
function at cu = also for interacting electrons;^ 

g N = ^4r i? (0)G+ 1 (0)r L (0)G^(0). (14) 



Here T L (u) = - Im \y\ gt{u)] , and T R (u) = - Im [v R g R {co)] . g\ (g%) is the Green's 
function at the interface i — (i — N + 1) of the isolated lead, and it is determined by 
Hl (Hr)- Note that Eq. fll4]) is also expressed in terms of the transmission probability 
of the free quasi-particles, \t(0)\ 2 = 4r fl (O)G+ 1 (O)r L (O)G'^ v (O).0iil Correspondingly, the 
reflection probability is given by |F(0)| 2 = 1 1 - 2zr L (0)G+(0)| 2 = |l - 2zr /? (0)G+ 7V (0)| 2 , 
and the unitarity |t(0)| 2 + |r(0)| 2 = 1 is preserved owing to the property Eq. ( |IUD [see also 
(p2])j. Another quantity which can be related to the scattering coefficients is the displacement 
of the total chargeillzl 



(15) 



Here (• • -) z and (••-)# denote the ground-state average of isolated leads described by Hl 
and 7Yr, respectively. At T = 0, AN tot can be expressed in terms of the S'-matrix for the 
quasi-particles, following the derivation of the Friedel sum rule by Langer and Ambegaokar,il 



as 



AN tot = — log[detS" 

7U 



(16) 
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(17) 



r L (o) o 
o r R (o) 

Therefore, the conductance and charge displacement are determined by the inter- and intra- 
boundary elements of the Green's function such as G^- 1 (0) and Gf^O). 

Next we discuss the structure of the Dyson equation Eq. further in order to make 
the mathematical features of the perturbation theory in the presence of the reservoirs clear. 
For the Green's functions in the central region, Eq. (P) is written in a N x N matrix form; 



{Q(z)}- 1 = {g®(z)}- 1 -V(z). 



(18) 



Here Q(z) = {Gij(z)} with ij £ C, and the inverse matrix of the unperturbed part can be 
expressed as {G^\z)}~ 1 = z 1 — H? c — V m i X (z) with 



L ll 



l 12 



n/0 



L 21 



L 22 



(19) 



"V mix \%) 

























E(s) 
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Sn(^) £i 2 (V) 
S 2 i(-2) £ 22 U) 







(20) 



(21) 



Here 1 is the N x N unit matrix. V m ix corresponds to the contribution of the mixing 
with the reservoirs. Especially, the two non-zero elements v\ g L and v 2 R g R have the finite 
imaginary parts corresponding to and T R defined just bellow Eq . ([HQ . These imaginary 
parts change the discrete levels of HP C to the continuous peaks with finite level width, and 
thus the unperturbed part of the Green's function Q^°\z) describes a system of iV resonant 
scatterers. Therefore, owing to the contribution of the mixing, the mathematical structure 
of the perturbation theory in the presence of the reservoirs becomes similar to that of the 
Anderson model with a number of orbits rather than the usual Hubbard model without 
reservoirs. In this matrix form, the quasi-particle description is summarized as follows. 
Due to the property ImS + (0) = at T = 0, the Green's function at u = is written as 
(£? + (0)} 1 = TC — V mix (0), where K is the renormalized hopping matrix corresponding to 
-He defined by Eq. (U); 



K = - [?4 + Re£ + (0)] . 



(22) 



Then a relation corresponding to the optical theorem can be obtained as G + (0) — Q (0) = 
£/ + (0) [V^ i;c (0) — V~ ix (0)] Q~(0), which certifies the unitarity of the scattering coefficients. 
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Here the superscript + (— ) means the retarded (advanced) function. 

The quasi-particle approach cannot be applied to a non-Fermi-liquid ground state where 
the perturbation theory with respect to the inter-electron interaction breaks down. However, 
in the finite interacting system connected to reservoirs as illustrated in Fig. [I], a cross over 
from a high temperature phase to the low temperature Fermi liquid phase is expected in 
many cases. This is because, due to the connection with the reservoirs, not only the energy 
scale corresponding to the level spacing of the isolated interacting system of size N but also 
the level width introduced through Tl and Tr, i.e., the imaginary part of the mixing term 
Eq. fl20D , plays a role of the cut-off. In this sense, as mentioned above, the system can be 
regarded as a generalized Anderson impurity with iV resonant states. 

When the interacting region is described by an one-dimensional chain, the Tomonaga- 
Luttinger (TL) behaviorS seems to be seen for large N. However, as discussed by Kane and 
Fisher for the finite TL model of size L connected to Fermi liquid reservoirs,! the cross over 
from the high temperature TL phase to the low temperature Fermi liquid phase occurs at a 
characteristic energy T ~ vp/L, where vp is the Fermi velocity. Since this simple discussion 
does not take into account various effects such as the back and Umklapp scatterings, the 
characteristic energy will be different depending on the details of the situations. Specifically, 
in the electron-hole symmetric case of the Hubbard chain examined in the following chapters, 
the characteristic energy depends on whether N is even or odd, and it should be the Mott- 
Hubbard gap E G or the Kondo temperature T K . Again, due to the contribution of the 
mixing, the low energy excitations below the characteristic energy scale are described by 
the Fermi liquid as far as N is finite. For even N, the Mott-Hubbard gap opens in the 
thermodynamic limit of the usual Hubbard model without reservoirs, and the corresponding 
tendency should be seen in the present system of finite N. However, there still remains 
a finite density of states at the Fermi energy caused by the mixing with the reservoirs, 
although the spectral weight in the region corresponding to the gap will decrease with 
increasing N. For odd N, there is an additional spectral weight at the Fermi energy in 
the Mott-Hubbard gap, i.e., the Kondo resonance of the width Tk- Thus, the insulator-like 
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low- conductivity should be seen at the temperature range Tk < T < Eg, and the metallic 
behavior of the Fermi liquid will be seen at low temperatures T < Tk- Note that Tk can 
be defined only for the open system connected to reservoirs, while Eq is a constant defined 
in the thermodynamic limit of the isolated system. Furthermore, Tk should decrease with 
increasing N, and in most of cases Tk should vanish in the limit of N — > oo. Therefore, 
if the limit iV — > oo is taken keeping T to be finite, the even-odd behavior disappears 
and at T < Eq the insulating behavior meets with an intuitive picture will be seen. This 
limit corresponds to the thermodynamic limit of the macroscopic system. However, there 
is another nontrivial limit which describes the low temperature physics of the mesoscopic 
systems, i.e., the limit T — > keeping iV to be finite. The ground state of the small system 
belongs to this limit, and in the case of the small Hubbard chain the even-odd behavior can 
be seen at low temperatures since Tk is finite for small N. This kind of a mesoscopic limit 
is expected to be realized in quantum dots or quantum wire of nanometer size, and our aim 
here is to develop a Fermi-liquid theory for the mesoscopic systems. 

Another possible non-Fermi-liquid state is the ground state of the multi-channel Kondo 
system.il Shimizu, Sakai and Suzuki have shown with the numerical renormalization group 
approach that the non-Fermi-liquid ground state is realized also in an extended version 
of the Anderson model.0 So, if the parameters of the Hamiltonian Eq. ([!]) satisfy certain 
conditions, the similar situation seems to be realized. 

III. QUASI-PARTICLES IN THE ELECTRON-HOLE SYMMETRIC CASE 

The effective Hamiltonian has some notable properties in the electron-hole symmetric 
case where the average number of electrons in each site is unity. In this section, we provide 
a simplified expression of the conductance in the electron-hole symmetric case. The results 
Eqs. (|31| ) and ([!3|) will be used for the finite Hubbard chain in the next section. Especially, 
form Eq. (|3T|), we can deduce quite generally that the perfect transmission occurs for odd 
N when the system has an additional inversion symmetry. 
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In the electron-hole symmetric case, the off-diagonal element of the renormalized pa- 
rameter tfj is zero when i and j belong to the same sublattice, i.e., —tf^ — fid^ = for 
N ~~ j\ = 0) 2, 4, . . ., and Re V^^O) = 0. Thus, the matrix 7C defined by Eq. ( p2|) has a 
checkered structure, and the Green's function at T = 0, to = is written in the form; 








7C 
' 12 





1 14 










7c 

L 21 





7C 
L 23 


■■• 




■ *r L (o) 




{e + (o)r= 





7c 
'32 







+ 





(23) 




t 41 





fc 43 


■■• 









and Gfj(0) is obtained by taking the inverse of Eq. fl23|). Specifically, using the explicit form 
of Gf N (0), the dc conductance g N can be obtained from Eq. (|TJ]). For N — 1 and 2, the dc 
conductance is written in the form 



2e 2 
9i — ~~r~ 



L r r 



2e 2 

92 = ~r 



h [(T L + T R )/2} 2 ' 

r L r R {t c 12 y 



(24) 



(25) 



h [(T L T R + { t C 2 y)/2} 2 
Here T L = Tl(0), T r = r^(0), and we will use this simplified notation in what follows. For 

N > 3, G± N (0) is written in the form 

det TC/vi 



\N+1 



det{S + (0)} 



-l ' 



(26) 



det{£ + (0)}^ 



(27) 



i T L det /On + i T R det K,nn , for odd iV 
-T L T R det TCf^ + det 7C , for even JV 
Here /Qj is a (N — 1) x (N — 1) matrix obtained from K. by deleting the i-th row and the 
j-th column. Similarly, /C^ is a (iV — 2) x (JV — 2) matrix obtained from 7C by deleting the 
first and the iV-th rows, and the first and the N-th columns. Due to the checkered structure 
of Eq. {\2,'S\), the determinants satisfy a relation 

det K n det/CiVAf , for odd N 



(det K N1 ) 2 



det 7C det K^ N , for even N 



(28) 
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Although the details are given in the appendix [A], the origin of the even-odd dependence 
can be understood just seeing typical examples; 



" 


t C 
L 12 





'11 
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TC 
L 12 





TC 

t 23 





TC 
t 52 





TC 
L 23 





TC 
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TC 
t-41 





rc 

'31 





TC 
''45 
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''45 






TC 





" 


TC 

L 12 





'11 





^61 




tl2C 





rc 

t 23 





TC 
L 52 










TC 
fc 23 





TC 
t 34 





TC 




TC 





TC 
t 34 





rc 

''45 










7C 
t 52 





TC 
1 45 





rc 




TC 





TC 
b 63 





TC 
1 56 





elements, 




-1-3 


for j 


= 1, 


2, .. 



for = 6 . 



(30) 



N, are zero for odd N, while those 
are not for even N. Furthermore, det TC is zero for odd N, while it is finite and can be 



divided into two determinants for even N. Using the properties Eqs. Q2TS|) 



(0), the dc conductance for odd N (= 2M + 1) is written in the form 

2e 2 f L f R 



#2M+1 



where T L = XT L , T R = T R / X, and 



A 



[(r L + r R )/2[ 



with Eq. 



(31) 



det/C 



li 



(32) 



det K-mn 

We note that the parameter A can be simplified as Eq. ( |A10| ) using the checkered structure 
further. If the system also has the inversion symmetry I\ = T R m addition to the electron- 
hole symmetry, the parameter is fixed as A = 1. This is because the matrix 7C becomes 
symmetric with respect to the antidiagonal line. Therefore, due to the combination of the 
symmetries, the perfect transmission occurs <? 2 m+i = 2e 2 //i independent of M and the details 
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of the interaction, as far as the perturbation expansion is valid. Physically, this is caused 
by the Kondo resonance appearing at the Fermi energy u — 0. On the other hand, for even 
iV (= 2M), the dc conductance is written in the form 

_ 2^ V L TrV 2 C 

ti [(r L r R + v£) /2\ 



vc = -^T^nn ■ (34) 



det 7C 
det/cj 7 

We note that the parameter vc can also be simplified as Eq. (|A6|) . These two expressions 
for the dc conductance, Eqs. (|3ll) and (|3~3"D, are the main results of this section. 

We next examine the Friedel sum rule. Using the property Eq. (]27|) and the expression 
of the total charge displacement Eqs. (|16|) and (|I7| ), we obtain det S = 1 for even N, 
and det S — — 1 for odd N. This is consistent with the fact that the average number of 
electrons in each site is unity. In order to realize the electron-hole symmetry, the bare matrix 
element tf, should also have the checkered structure as that is assumed for the renormalized 
parameter tfj, i.e., the system must be classified into two sublattices. Note that, since only 
the symmetry is assumed so far, the bare hopping matrix element is not necessary to be 
restricted to the nearest-neighbor one. Furthermore, the system is possible to be disordered 
through the randomness in the off-diagonal element of tfj. 

IV. TRANSPORT THROUGH A FINITE HUBBARD CHAIN 

In this section, we will use the effective Hamiltonian to investigate the conductance of 
a finite Hubbard chain connected to two semi-infinite leads. The system is considered as 
a model for a series of quantum dots or atomic wire of nanometer size. We calculate the 
renormalized hopping matrix element tfj up to the second order in U, and obtain the dc 
conductance using the expressions given in the previous section. We note that the second- 
order perturbation theory has been used by several groups for studying transport properties 
of a single impurityiB~H and systems consisting of a number of resonant levels.B^ For 
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N > 2, the off-diagonal part of the self-energy plays an important role on the conductance 
at T = as seeing in the results presented below. To our knowledge, however, effects of the 
off-diagonal part have not been examined sufficiently so far. 

The schematic picture of the model is illustrated in Fig. |^: the system consists of N 
interacting sites at the center 1 < % < N, and two non-interacting leads at — oo < i < 
and A^+l<i<+oo. The explicit form of the Hamiltonian is given by 

H = H + H I} (35) 

+oo +oo 

#0 = - X) *< ( c *+^ c ^ + c ^ c *+^ ) ~ ^ Yl Yl c ^ c ™ 

i= — oc cr i= — oo a 

+ EEG°+fK> < 36 > 



i=l a 
N 



i=l 



1 



n i\ n n - 2 K + n u 



(37) 



where U is the repulsive interaction in the chain, eo is the on-site energy, and n if7 = c ia c ir7 . We 
take the hopping matrix element to bet i = t for every links except for the coupling between 
the chain and leads, i.e., t = v L and t N = v R . Thus, the level width caused by the coupling 
with the leads is given by r a (0) = 7ri&D(0) for a = L, R with D(Q) = ^/4t 2 - // 2 / (2irt 2 ) . 
In what follows, we concentrate on the electron-hole symmetric case taking the parameters 
to be \i = and e = —U/2. 

Within the second order, the self-energy correction is described by a diagram Fig. [|, and 
the value at T = 0, to = is obtained from the expression 

S<?(i0 + ) = -V* ££ Gf(ie) G<°>(«<) G<?( i£ + ie>) , (38) 

for 1 < i, j < N. Here G^\ie) is the unperturbed Green's function corresponding to 
H , and thus E^(i0 + ) depends on v L and v R through G^. Note that the value of the 

(9) 

retarded function at to = can be obtained from the Matsubara function, i.e., E^- (i0 + ) = 
S^(ie)| e _ > o+- The imaginary part of the self-energy vanishes owing to the Fermi-liquid 
property, i.e., ImE^ (z0 + ) = 0. Furthermore, n[f(i0 + ) = for even — since G^\—ie) = 
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(— Gtj(ie) in the electron-hole symmetric case. We have done the integration Eq. 
(j38[) numerically to estimate the renormalized matrix element tfj defined by Eq. (13) within 
the second order, and obtain the dc conductance g N using Eqs. (|3T]) and fl33|). As an 
example, we show the result of the self-energy for N = 6, 



2nt J 



(39) 



,(2) 





0.8596 


-0.1642 


0.0783 



0.8596 


0.3969 


-0.0352 






0.3969 


0.7543 


-0.1642 



-0.1642 


0.7543 


0.3969 






-0.0352 


0.3969 


0.8596 



0.0783 


-0.1642 


0.8596 




(40) 



Here the coupling is taken to be v L /t = v R /t = 0.8. The matrix element for a given 

(2) 

distance \i — j\ shows an oscillatory behavior, i.e., S ; i i f2m+1 (i0 + ) is an oscillatory function 



of I. In contrast, in the special case v. 



'R 



t, £ 



(2) 



'■> ^//+2m+i(^ + ) is independent of I and 
the integration Eq. (38) can be done analytically.!^ Note that T,[J ) (i0 + ) is also an oscillatory 
function of \i — j\, and the absolute value tends to decrease with increasing \i — j\. 



A. Symmetric connection vl = Vr 

We now examine the case where the system has the inversion symmetry, v L = v R = v, 
in addition to the electron-hole symmetry. In this case, as it was shown generally in the 
previous section, the perfect transmission occurs for odd N independent of U. Physically, 
this is due to the contribution of the Kondo resonance appearing at the Fermi energy. So, we 
have done the numerical calculation only for even N . The result of g N is plotted as a function 
of the size N in Fig. ^, where the parameters are taken to be v/t = 0.8 and U/(2nt) = 0.4. 
The result shows a typical even-odd oscillatory behavior, and the conductance for even N 
decreases with increasing N. This reduction is caused by the oscillatory I dependence of 
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the self-energy S}/ 2m+ v In Fig. |^, the envelope of g N for even A" (= 2M) is plotted vs M 
for several values of the repulsion U/(2wt) = 0.0, 0.2, 0.4, 0.6, and 0.8, taking the mixing 
matrix element to be v/t = 0.8. The conductance g 2M decreases with increasing M and 
increasing U. In the non-interacting case, g 2M is independent of M since the parameter 
defined by Eq. (|34] ) is simply vc = t. The reduction of the conductance seems to show a 
tendency toward a Mott-Hubbard insulator. However, the result shows a rather moderate 
decay, i.e., it does not show an exponential dependence. Since the range of U in which the 
second-order perturbation theory is able to provide quantitatively reliable results depends 
on the size M, contributions of higher-order terms should be examined in order to clarify 
the correct asymptotic behavior for large M. 

In Fig. ^, the conductance is plotted as a function of U for a number of even iV (= 
2, 4, 6, . . .) taking the parameters to be v/t — 0.8 (dashed lines) and v/t = 1.0 (solid lines). 
As it can be seen in the behavior of the dashed lines, the reduction of g 2M is proportional 
to U 2 for small U when \v/t\ < 1.0. The curvature increases with the size M. In order to 
see this quantitatively, we expand the conductance in powers of U as g 2M = (2e 2 /h)[C - 
C 2 {U / (27rt)} 2 + •••], and plot the ratio C 2 /C as a function of M in Fig. [7| for several values 
of v/t (= 0.1, 0.2, . . . , 0.9). C 2 increases with the size M, and seems to diverge in the limit 
of large M showing a power low behavior. Note that the transmission probability in the 
non-interacting case Cq depends on the mixing matrix element v but is independent of the 
size M as mentioned in the above. Furthermore, since the self-energy is calculated up to the 
order U 2 , the result of C 2 is exact. The ratio C 2 /Cq increases with decreasing v. This means 
that the effect of the interaction is enhanced when the coupling between the sample and 
leads is weak. In the special case v = t, the reduction of g 2M is proportional to U A for small 
U, i.e., g 2M = (2e 2 /h)[l — C±{U '/ '(27rt)} 4 + •••]. This is because the unperturbed system 
described by Hq has a translational invariance accidentally in this case and the reflection 
probability is zero at U = 0. The solid lines for > 10 in Fig. ^ are almost overlapping 
each other. In Fig. |8|, C4 is plotted vs M. This result is also exact, and the coefficient C4 
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converges to a finite value for large M; C4 ~ 0.5293 for N = 250. In the appendix |B|, the 
convergence of C4 in the limit of M — > 00 is confirmed with another approach. Therefore, 
when v = t, g 2M is finite even in the limit of large M at least for small [7.0 For large values 
of U, both the dashed and solid lines in Fig. |^ tends to zero showing a l/U 4 dependence. 
However this behavior depends on the approximation: as it can be confirmed from Eqs. 



and (34|), vc oc U n and g 2M oc 1/U 2n for large U when the self-energy is estimated within 
the n-th order perturbation in U. 

The mixing matrix element v determines the bare level width of the resonant states, 
and in the present case it is given by T = v 2 /t. Since we are now considering the electron- 
hole symmetric case, the Fermi level for even N is located between the two resonant states 
corresponding to the highest occupied and lowest unoccupied levels. In Fig. |9|, g 2M is plotted 
vs U for several values of v/t (= 0.4, 0.5, 0.6, . . . , 1.0), where the solid (dashed) lines are 
the results for N = 40 (N = 4). The value of the conductance itself decreases with v, but 
the qualitative feature of the U dependence is similar in the cases of \v/t\ < 1.0. In Fig. 
[TP] , we show the conductance as a function of M for several v/t (= 0.4, 0.5, . . . , 1.0) taking 
the repulsion to be U/(2irt) = 0.3. In this figure, the conductance is normalized by the 
non-interacting value g^M = (2e 2 /h) Cq which depends on v. The normalized conductance 
decreases with v. This also means that the effect of the interaction is enhanced when the 
level width of the resonant states is small. 



B. Asymmetric connection vl 7^ vr 

We next examine the case where the inversion symmetry is broken, v L 7^ v R , but the 
electron- hole symmetry is still preserved in the equilibrium state by the condition \x = 
and €q = —U/2. As one typical example, we take the mixing parameter to be v L /t = 0.8 
and v R /t = 0.6 in this subsection. In Fig. [XT] (a) and (b), the conductance is plotted as a 
function of the size N taking U/(2wt) to be (a) 0.2 and (b) 0.4, respectively. The result 
shows the even-odd oscillatory behavior as that in the case of v L = v R discussed in the 
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above. However, in the case of v L 7^ v R , the perfect transmission does not occur as it can 
be deduced from Eq. (|3~lD, so that we have done the numerical calculation also for odd 
N. Although the Kondo resonance is still present at the Fermi level owing to the electron- 
hole symmetry, the conductance is reduced. Nevertheless, the Kondo state contributes to 
the higher transmission for odd N. In Fig. [12], g N is plotted as a function of U for a 
number of iV (= 1, 2, 3, . . .,40). For even N, qualitative features of the results are similar 
to those in the inversion symmetric case [See Fig. On the other hand, for odd N, the 
features of the curves are somewhat different. However, the results for large values of U 
have some limitations, because the parameter A defined by Eq. (|32D tends to an incorrect 
limit for large U when the self-energy is estimated within a finite order in U, i.e., A tends 
to the ratio of the highest-order term of the numerator and that of the denominator of Eq. 
( |A10p . In the non-interacting case the conductance for even and odd N, g 2M and g 2 M+u 
are independent of the size M, respectively, since vc — t and A = 1. For small U, the 
reduction of the conductance is proportional to U 2 and the curvature increases with the size 
M, except for the single impurity case N = 1. In Fig. the coefficient C 2 defined by 
g N = (2e 2 /h)[C — C 2 {U / (2irt)} 2 + ■ ■ ■] is plotted as a function of N. In the figure, the 
envelop for even N is in upward from that for odd N. The coefficient C 2 increases with the 
size N showing an oscillatory behavior, and probably it diverges in the limit of large N. 

V. SUMMARY 

We have applied a quasi-particle description of a Fermi liquid to the transport through 
a small interacting system connected to reservoirs. With this approach, we have studied 
the properties of quasi-particles in an electron-hole symmetric case. In this case, effects 
of electron correlation come in the theory through the off-diagonal part of the self-energy, 
and the conductance can be written in a simplified form Eqs. (|3~1~D or (|33|) depending on 
whether the number of the interacting sites N is even or odd. It is shown using Eq. fl3T|) 
that the perfect transmission occurs quite generally for odd N if the system has an inversion 
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symmetry in addition to the electron-hole symmetry. We have applied the method to a 
small Hubbard chain of finite size N with the second-order perturbation theory in U, and 
examine the conductance as a function of N, U, and the coupling between the chain and 
the leads, i.e., v L and v R . Effects of electron correlation are enhanced in weak coupling case 
v L , v R < t. In this case, the conductance for even N tends to zero for large N. This behavior 
is qualitatively different from that in the special case v L = v R = t. This seems to be caused 
by the difference in the structure of the resonant states. Quantitatively, the range of U in 
which the second-order perturbation theory is able to provide reliable results tends to be 
narrow with increasing N, and thus the higher order terms should be included in order to 
refine the results obtained for large U and to clarify the asymptotic behavior for large N. 

Throughout this work, we have assumed that the interaction is switched on only in the 
sample region at the center. When the interaction is switched on also in the two leads, 
the description in terms of the quasi-particles is still possible if the ground state is changed 
continuously against this interaction. But some modifications are necessary. The interaction 
in the leads will cause the renormalization of the incoming and outgoing fields corresponding 
to initial and final states of the scattering matrix, i.e., the renormalization of the external 
lines. Furthermore, the quasi-particle description can be extended to finite temperatures by 
taking into account the residual interaction among the quasi-particles or the contributions 
of vertex corrections. Especially, an overall picture of the even-odd property described in 
the last part of Sec. |T| will be confirmed microscopically based on the finite temperature 
theory. 
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APPENDIX A: 7C MATRIX IN THE ELECTRON-HOLE SYMMETRIC CASE 



We summarize here properties of the renormalized hopping matrix 7C, and provide the 
derivation of Eq. fl5Sp. In the electron-hole symmetric case, K has a checkered structure 
as shown in Eqs. fl29|) and fl30|). Thus, det K = for odd N. For even N (= 2M), the 
determinant can be divided into two parts; 






rc 

L 12 





rc 
L u 


fC 
'21 





rc 

'23 








rc 





rc 

^34 


''11 





rc 

' 13 






;-l) M detQ det Q' 



(Al) 



Here Q and Q' are M x M matrices defined by 

r 7c 



12 



rc 

'32 



rc 

''31 



fc 21 
t 41 



^23 
fc 43 



(A2) 



Since the renormalized parameter is real and tfj = owing to the time-reversal symmetry, 
det Q = det Q' . Consequently, for even A" (= 2M), the determinant of the matrices K, 



K U N , and Km can be factorized as 



det K 

•NN 



-l) M {detQ} 2 , 
det 70^ = (-l)^- 1 {detQf^} 2 
det K N1 = det Q det Q[\ N . 



(A3) 
(A4) 
(A5) 



Here Q[\ N is a (M — 1) x (M — 1) matrix extracted from in the similar way as it was 
done for extracting Q from K. Using Eqs. (|A3|) -([A5|), we obtain Eq. (p8|) for even N, i.e., 



(det/C 



Nil 



— det K det K xx . Furthermore, vq defined by Eq. ([34] ) is simplified as 



vc 



detQ 
detQ™ 



(A6) 



For odd A" (= 2M + 1), the 2M x 2M matrices TCn and Knn have the similar properties 
described in the above. Consequently, the determinants can be factorized as 
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det /Cn 



{detQ n } 2 



detK NN =(-l) M {det Q NN Y , 
det Kni = det Q n det Q NN . 



(A7) 
(A8) 
(A9) 



Here Q n and Q NN are M xM matrices extracted from TCn and JCnn, respectively, as it was 
done in the above. Using Eq. (|A^-(|A|), we obtain Eq. @ for odd N, i.e., (detK N1 ) 2 = 
det /C 11 def/C^vv- Also, A defined by Eq. (|32|) is simplified as 

det Q u 



A 



detQ 



NN 



(A10) 



APPENDIX B: REFLECTION COEFFICIENT IN A SPECIAL CASE 

We show here another approach to clarify the asymptotic behavior of the coefficient C4 
introduced in Sec. [TV] in the case of v L = v R = t. In this case, the Dyson equation Eq. (g) 
is written in terms of the scattering matrix Tu> as, 

N 

G,r = G«°! + V Gf %, Gf], . (Bl) 
M'=i 

The transmission and reflection coefficients are also written in terms of Tw as t(0) — 1 — 
iT kFkF {iO+)/v F and r(0) = -i T_ kFkF {tO + )/v F £M where v F = 2t, k F = vr/2, and 

N 

T kk ,= e~ lkl T w e ik ' 1 ' . (B2) 
i,v=i 

Since the unperturbed system has the translational invariance in the present case, the lowest- 
order scattering matrix is given by T^p = E^, . Thus, using a wavenumber representation 
of E^S the lattice sum Eq. ( |B2| ) can be done explicitly, and the reflection coefficient is 
expressed as 

dk s dk 2 dk 1 cos k F N — cos[(fci + k 2 — k^N] 



^2) L _ e ik F {N+l) jj2 



A v F J (27r) 3 cos k F — cos {k\ + k 2 — k 



3 J 



x p /fc 3 (l - /fc 2 )(l - /fci) + (1 - fk 3 )fk 2 fk 1 ^gg^ 
£,k 3 S&2 S&l 
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Here = — 2tcosk, = [e^ k + and P denotes the Cauchy principal value. It can be 
confirmed from Eq. ( |B3| ) that = for odd N. For large even N, the contribution of the 
fast varying cos[(fci + k 2 — ks)N] part becomes small, and the dominant contribution comes 
form the cos kpN part. Thus, for large even N (= 2M), the reflection coefficient is written 



in the form lim r% M = V^T {U/^nt)} 2 with 

M— >-oo 

4 An y_ 7r / 2 [cosfci +cosA;2 H-cos/^] cos (k\ + k 2 + k 3 ) 



The coefficient -^/Cf is finite and estimated numerically as \[C^ — 0.729. Thus, in the 
limit of large even N, the dc conductance in the case of v L — v R — t is written in the form 
lim g 2M = (2e 2 /h) [ 1 - C 4 °° {U/(2itt)Y + ■ • • ] . 
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FIGURES 



FIG. 1. Schematic picture of the system 
FIG. 2. Schematic picture of the model: (•) interacting region, (o) ideal leads. 

(2) 

FIG. 3. Second-order self-energy EJ- (iu). 

FIG. 4. Conductance g N as a function of the size N in the inversion symmetric case. Here 
v/t = 0.8 and U/(2wt) = 0.4. 

FIG. 5. Conductance for even N (= 2M) as a function of the size M. Here v/t = 0.8, and the 
repulsion U/(2irt) is taken to be (•) 0.0, (■) 0.2, (♦) 0.4, (▲) 0.6, and (*) 0.8. 

FIG. 6. Conductance vs U, for even N (= 2, 4, 6, . . .). Here v/t = 0.8 for dashed lines, and 
v/t = 1.0 for solid lines. Totally, 20 dashed and 30 solid lines are plotted. 

FIG. 7. The ratio C2/C0 is plotted for several v/t (=0.1, 0.2, . . . , 0.9) as a function of the size 
M. Here g 2M = (2 e 2 /h)[C - C 2 {[//(2vrf)} 2 + •••]. 

FIG. 8. The coefficient C4 in the case of v/t = 1.0 is plotted as a function of the size M. Here 
ff2M = (2e 2 // l )[l-C 4 {[//(2vrt)} 2 + ...]. 

FIG. 9. Conductance vs U, for N = 4 (dashed lines) and N = 40 (solid lines). Here v/t is 
taken to be 0.4, 0.5, . . . , 0.9, and 1.0. 

FIG. 10. Conductance for even N (= 2M) as a function of the size M. Here U / {2irt) = 0.3, 
and the conductance is normalized by the value for U = 0, i.e., = (2e 2 M)Co. ^ ne mixing 
matrix element v/t is taken to be (■) 0.4, (•) 0.5 (*) 0.6, (▲) 0.7, (♦) 0.8, (■) 0.9, and (•) 1.0. 

FIG. 11. Conductance g N as a function of the size N in the inversion asymmetric case. Here 
v Jt = 0.8, v R /t = 0.6, and U/{2-Kt) is taken to be (a) 0.2 and (b) 0.4. 
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FIG. 12. Conductance vs U, for a number of N (= 1, 2, 3, ...,40). Here vjt = 0.8 and 
v R /t = 0.6. 

FIG. 13. The coefficient C2 as a function of the size N. Here v L /t = 0.8, f^/t = 0.6, and 
g N = (2e 2 /h)[C - C 2 {U / (2irt)} 2 + ■••]. 
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